Image Properties of Embedded Lenses 

R. Kantowski, 1 *!] B. Chen/fl and X. Dai 1 ^ 

1 Homer L. Dodge Department of Physics and Astronomy, 
University of Oklahoma, 440 West Brooks, Norman, OK 73019, USA 

(Dated: March 8, 2013) 



1 



Abstract 

We give analytic expressions for image properties of objects seen around point mass lenses em- 
bedded in a flat ACDM universe. An embedded lens in an otherwise homogeneous universe offers 
a more realistic representation of the lens's gravity field and its associated deflection properties 
than does the conventional linear superposition theory. Embedding reduces the range of the grav- 
itational force acting on passing light beams thus altering all quantities such as deflection angles, 
amplifications, shears and Einstein ring sizes. Embedding also exhibits the explicit effect of the 
cosmological constant on these same lensing quantities. In this paper we present these new results 
and demonstrate how they can be used. The effects of embedding on image properties, although 
small i.e., usually less than a fraction of a percent, have a more pronounced effect on image dis- 
tortions in weak lensing where the effects can be larger than 10%. Embedding also introduces a 
negative surface mass density for both weak and strong lensing, a quantity altogether absent in 
conventional Schwarzschild lensing. In strong lensing we find only one additional quantity, the 
potential part of the time delay, which differs from conventional lensing by as much as 4%, in 
agreement with our previous numerical estimates. 
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I. INTRODUCTION 

Recently we have investigated the quantitative effect of embedding on gravitational lens- 
ingobservations by resorting to a mixture of analytic work with a few numerical applications 
[l-3]. The analytic results for quantities like the bending angle a produced by a point mass 
were given as functions of two impact variables r and <j)\ (see Fig.l). These two parameters 
are not independent if the source and deflector redshifts are fixed. Because of the non- 
linearity of the expressions we were only able to give an iterative procedure that allowed us 
to numerically evaluate the conventional minimum impact Schwarzschild coordinate tq as a 
function of 0i We have since been able to analytically carry out this iterative procedure 
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see Eq. ( lA.lj) in the appendix) and hence obtain all lensing properties such as position, 



shear, etc., as functions of the single impact angle <j>i. The solution of the embedded lens 
equation and comparison with classical lensing theory is therefore greatly simplified. Be- 
cause the dependence of observable quantities on this angle is highly nonlinear, we are not 
able to eliminate <f>i in favor of r . Derivations of our current results follow the steps given 



in 



u 



3| which we will not repeat but we will instead simply present the new results and use 
them on two examples. 

The point mass lens is the simplest lens to use to demonstrate the effects of embedding; 
however, all lenses will require corrections. An embedded point mass lens is constructed by 
condensing a comoving sphere of pressureless dust of a standard homogeneous cosmology to 
a singular point mass m at the sphere's center, a construction first made by Einstein himself 



0£ 



HJ|. When the ecology contains a cosmological constant A the gravity field inside the 
evacuated sphere is described by the Kottler metric [8|, [9( rather than the Schwarzschild 
metric. In this paper we restrict ourselves to a flat background cosmology whose Friedman- 
Lemaitre-Robertson- Walker (FLRW) metric is 

ds 2 = - c 2 dT 2 + R(T) 2 [d X 2 + X 2 (d8 2 + sin 2 6d<j) 2 )] . (1) 
The embedded condensation is described by the Kottler or Schwarzschild-de Sitter metric 
ds 2 = - 1 {r)~ 2 c 2 dt 2 + 7 (r) 2 rir 2 + r 2 (d9 2 + sin 2 9 d<j) 2 ), (2) 

where 7~ 1 (V) = y/l — fi 2 {r) and /3 2 (r) = r s /r + Ar 2 /3. The constants r s and A are the 
Schwarzschild radius (2Gm/c 2 ) of the condensed mass and the cosmological constant re- 
spectively. By matching the first fundamental forms at the Kottler-FLRW boundary, angles 
(9, 4>) of equations (JT|) and ()2]) are identified and the expanding Kottler radius of the void 
is related to its comoving FLRW boundary Xb by 

r b = R(T) Xb . (3) 

By matching the second fundamental forms the comoving FLRW radius Xb is related to the 
Schwarzschild radius r s of the Kottler condensation by 

r s = n m ^(R oXb ) 3 . (4) 

Here Hq is the familiar Hubble constant and the cosmological constant A is constrained to 
be the same inside and outside of the Kottler hole. 
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FIG. 1. A photon travels left to right entering a Kottler hole at r = r%, (f> = n — <p\ and returns to 
the FLRW dust at r = r\ + Ar, <j) = (f>i + A<^>. The photon's orbit has been chosen symmetric in 
Kottler about the point of closest approach r = ro, 4> = ir/2. Due to the cosmological expansion, 
Ar > 0. The slope of the photon's co-moving trajectory in the x-y plane is £i when incoming 
and £i + a after exiting. The resulting deflection angle as seen by a comoving observer in the 
FLRW background is a, which is negative by convention. Expressions for r±, Ar, £i, A</>, and a as 



functions of the two impact parameters, ro and (f>i, can be found in 
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In the following sections we will give image locations and image properties of small sources 
seen through Kottler voids in an otherwise flat FLRW universe (an embedded lens). We 
assume that the source and deflector are located at fixed FLRW comoving distances Xs 
and Xd from the observer which correspond to angular diameter distances D s and Dd, and 
redshifts l + z s = Rq/R s and 1 + Zd = Ro/Rd, see Fig. 2. These quantities are computed just 
as if the void didn't exist. Any quantity with a subscript '<f means that it is evaluated at 
redshift Zd when the radius of the universe was Rd = R(Td) = Rq/(1 + Zd). We give lensing 
properties such as the bending angle a of Eq. (JZJ) that are a series of smaller and smaller 
terms, sufficient to see both the shielding effect of embedding and the effect of the expansion 
rate (3d = Vd/c of the void's Kottler radius = RdXb that existed at FLRW time Ta- The 
expansion rate Vd is the speed of the expanding void boundary as measured by a stationary 
Kottler observer at r^. It is given by evaluating (3(r) defined below Eq. ([2]) at r — 



d 




(5) 



When expanding quantities such as a in a series we have taken parameters (3d and Xb/Xd 
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FIG. 2. A photon travels from a source at comoving distance Xs from the observer and then enters 
a Kottler hole of comoving radius Xb centered at comoving distance Xd from the observer. The 
photon is deflected by an angle a < and returns to the FLRW dust on its way to the observer. 
Because this is a comoving picture the orbit inside the void is only representative and because 
the true orbit inside the void is symmetric about = tt/2 (see Fig. 1) the optical axis is rotated 
clockwise by the angle — p (see Eq. (14) of j]]). 

TABLE I. Embedded lens parameters for two deflectors (a galaxy and a cluster) at z d = 0.5 when 
viewing a source at z s = 1.0 in a J7 m = 0.3, Qa = 0.7 universe with Hq = 70kms~ 1 Mpc -1 . 
Choosing R = 1, gives Xd = 1-89 x 10 3 Mpc and Xs = 3.31 x 10 3 Mpc. 

Lens m [3 d r d /D d = Xb/Xd r s /r d Ar^/3 ^(rad) </>£(rad) 

galaxy 10 12 M Q 3.67 x 10~ 4 9.54 x 10" 4 7.98 x 10 -8 5.51 x 10 -8 8.07 x 10" 6 8.45 x 10 -3 
cluster 10 15 M(T) 3.67 x 10 -3 9.54 x 10 -3 7.98 x 10 -6 5.51 x 10~ 6 2.55 x 10" 4 2.65 x 10~ 2 



r d /D d (the angular radius of the Kottler hole, see Fig. |2]) to be first order and r s /r d and 
Ar^/3 to be second order. In our results, e.g., Eq. ([7]), we have used a parameter 5 to keep 
track of each order. In Table 1 we give values for these and other parameters for two lens 
masses, m = 1O 12 M (a large galaxy) and m = 1O 15 M (a rich cluster) both at redshift 
z d = 0.5 in a flat Q m = 0.3, = 0.7 universe with Hq = 70kms _1 Mpc -1 . We refer to 
these as the galaxy lens and the cluster lens throughout the paper. 

Shielding typically causes the most significant embedding effect on images (i.e., the lowest 
order effect) and analytically appears as combinations of trig functions of the impact angle 
<pi in quantities like the bending angle a (see the first cos 3 0i term in Eq. (J7j) below). This 
decrease in a is caused by the shortened period a passing photon is influenced by the 
mass condensation (recall that in conventional lensing the deflecting force has "oo" range). 
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Corrections caused by the presence of A first appear in the void's expansion rate (3d and 
are typically smaller than shielding; corrections (i.e., are higher order). It was the search 



or A's effect on light deflections 
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16| that prompted investigations of embedded lensing 



II. THE EMBEDDED LENS EQUATION 

In our results we have introduced an order parameter 6 whose value is equal to 1 but 
whose purpose is to keep track of terms of similar orders as defined in the previous section. 
For weak lensing impact angles 0i, the higher the power of 8 the smaller the respective 
terms. For strong lensing, when 0i is sufficiently small, not all terms of a given order are 
of the same magnitude. By using steps developed in jl-^] we find that to order 6 4 the 
source and image positions 8s and 9j, as functions of the single impact parameter 0i, can 
be written as 

Xb 



g s = e I + ^a 



2(Xs - Xd) 



— sin 2 0i H — (3d I 4 — sin 2 0i + 3 sec 0i lo£ 

Xd 3 



(see Fig. 2), where the bending angle is 
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We identify the pair of equations <Q and ([8]) above as the embedded lens equation in 
parametric form. They can be used to obtain image positions and properties just as in 
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FIG. 3. Primary and secondary image positions as functions of source position 8s for the cluster 
lens of Table 1. Beyond 6 s ~ 3A6e the smallness of the secondary image's impact begins to violate 
the orbit approximation condition sin</>i 3> r s /rQ. 

conventional lensing theory. The conventional non-embedded lens equation is recovered 
by keeping only the lowest order terms in each expression and assuming cos^i —> 1 and 
r d sin0i r . 

Because of the dependence of each 8 order on the impact angle 0i, higher order terms 
that contain trig functions like esc <pi or cot <pi can become comparable in magnitude with 
the next lower order terms for sufficiently small values of <f)\. This happens in strong lensing. 
For example the embedded Einstein ring size 9' E is found by first finding the value of <pi = 4>e 
that makes 9s of Eq.((6]) vanish (see <j)E values in Table 1) and then evaluating 9' E = 9i(4>e) 
using Eq. ([8]). For 9s to vanish, S and <5 2 terms must cancel. The result is 

{9' E f = 9% cos 4>% (l + 38/3 d tan 4> E sin 4> E + C(<5 2 )) , (9) 

where 9e is the conventional Einstein ring radius defined by 

el = ^f. (io) 

As we have found with most strongly lensed image properties, this value differs only slightly 




from the conventional value. For the Einstein ring radius the embedded value differs some- 
what more than 0.05% for the cluster lens and 0.005% for the galaxy lens. In Figure 3 
we have used the embedded lens equation to locate primary and secondary images for the 
cluster lens. Primary and secondary images positions are given by Eq. (JSJ) and correspond 
respectively to impact angles <p + and 0_ (the Einstein impact angle <pE separates the two 
image domains, i.e., 0_ < 4>e < 4>+)- For a given source position 0$, primary and secondary 
image impact angels (fi± are found by solving 6s(<p±) = ±Os (i.e., by inverting Eq. (JSJ)). The 
two images are then located at ±0/(0±) (i.e., by using Eq. (jSJ). These two values of 0i can 
then be used to determine primary and secondary image properties. 



III. IMAGE PROPERTIES OF THE EMBEDDED LENS 



To evaluate standard image properties the reader only has to compute the azimuthal and 
radial eigenvalues (a^, a r ) of the image matrix dO s /dOj using equations (jSJ) and (jHJ). We give 
them in Equations ( 1 A. 21) and (1A.3I) of the appendix. The primary and secondary values for 
dcf, and a r can then be used to obtain image am plif ication, effective surface density, shear, 
and eccentricity, respectively fi, k, 7 s , and e (see 2(], 21]) by evaluating 




(11) 
(12) 
(13) 

(14) 

The above expressions give image properties for all values of impact angle <p\ such that 
the photon's orbit approximation is valid (sin 2 0! ^> r s /r^), but because of the lengths of 
the resulting expressions, we find it appropriate to make two approximations in the next 
section, one for weak lensing and one for strong. The effective surface mass density k for 
the embedded lens is the one property that does not vanish as it does for the conventional 
Schwarzschild lens and is plotted in Figure 4 for both weak and strong lensing of the primary 
cluster image. By a conventional Schwarzschild lens we mean conventional linear lensing 
theory applied to a point mass superimposed on a FLRW background. Even for strong 
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FIG. 4. The effective surface mass density k for the primary image of the embedded cluster lens 
of Table 1. 

lensing by the cluster the magnitude of k is only ~ 0.1% of the critical value. For the galaxy 
lens a k plot is similar to the cluster plot but is approximately a factor of 10 smaller. 

IV. WEAK AND STRONG APPROXIMATIONS 

We found it necessary to keep terms to order d 4 in expressions such as a and 9j to 
obtain sufficiently accurate results for most strong lensing quantities. Most weak observable 
quantities do not require such accuracy. By dividing the domain for <f>i into strong and weak 
parts we are able to give shorter expressions for the two eigenvalues (a^, a r ) of equations 
flA.2j) and flA.3j) and hence simpler expressions for /i, etc. For the strong domain we take 
O.40£ < 0i < 5<pE and for the weak 50£ < 0i < tt/2. The maximum value for 6j is 
approximately the ratio r^jDd which from Table 1 is ~117 times the Einstein ring radius @e 
for the cluster lens and ~37 times for the galaxy. Strong lensing consequently occurs for Qj 
values up to ~5#e and weak lensing begins to occur when 6j exceeds that value. To obtain 
shortened expressions for weak lensing we need only keep terms of order S 2 . This allows 
us to determine the lowest order effects of lens shielding and void expansion (the /3d term) 
on image properties in the weak domain 5(/>e < 4>i < tt/2. We find that the approximate 
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expressions are accurate to at least 0.1% down to 0i = 5 0g for the cluster and to at least 
0.03% for the galaxy. For weak lensing Eqs. (1A.2|) and ( 1A.3|) simplify to 

a ^ e afc = 1 _ 6 (2e£±^ csc 2 0icos 3 1 |l + (5/3 ( isec0i(2 + sin 2 1 ) + C(5 2 )|, (15) 
a weak = 1 + g (^E^j esc 2 0i I cos 0i (1 + 2sin 2 0i) + 8f3 d {2 - sin 2 0! + 2sin 4 0i) + C(<5 2 )|. 
From these the following image properties result: 

^eakyl = 1+3 $ (^Rl^j I COS 01+ 6/3 d sin 2 01 + O (S 2 ) 

-S 2 4 cot 4 0i I (1 + 2 sin 2 0i) + O (5) | , 

(^r) 2 { cos ^ + s & sin2 ^+o{8 2 )^ 

S (^^j esc 2 0i|cos0i + ~ sin 2 0^ + 6f3 d (2 - sin 2 0i + i sin 4 0^ + (<5 2 ) 

e" eafc = jv 7 ^ ) csc0iy / 2cos0i(2 + sin 2 0i) + £ /3 d (4 - 2 sin 2 0i + sin 4 0i) + 0(<5 2 ) 

xl/jl + <5^^ esc 2 0i cos 0i (1 + 2sin 2 0i) + 0(S) 

In Figure 5 we have compared image shear and ellipticity of the embedded lens images 
with conventional (non-embedded) Schwarzschild values. The reader can see that beyond 
0i ~ 45° the embedded lens differs from Schwarzschild by over 10%. This is caused pri- 
marily by the shielding of the embedded mass and increases as the transiting light ray's 
minimum impact r approaches the void boundary. The embedded amplification /1 differs 
from conventional Schwarzschild by less than 0.2% for the cluster and 0.02% for the galaxy 
for the weak lensing domain and only increases to 2.5% for the secondary cluster image 
in the strong lensing limit 0i — > 0.40^ where fi — > —0.03. The galaxy lens' numbers are 
significantly less and neither are plotted. 

Using the order parameter 6 to track terms of equal importance is problematic for strong 
lensing. As discussed above for strongly lensed images, small values of 0i in trig functions 
like esc 2 0i increase the numerical magnitudes of some of the terms in Eqs. (1A.2|) and (1A.3I) . 



In the following strong lensing approximation we have kept terms based on their numerical 
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FIG. 5. Corrections to image shear (left panel) and ellipticity (right panel) caused by embedding 
of the cluster lens of Table 1. The solid blue curves in the left and right panels are respectively the 
7s and e for the conventional (non-embedded) Schwarzschild lens. The fractional difference in the 
image shear, A7 S /7 S , and ellipticity, Ae/e, caused by embedding, are the dashed red curves in the 
left and right panels, plotted as a functions of <fr\. Differences are computed at the same primary 
image positions #r(0i). 



size at the Einstein ring value 0i = 4>e, and ordered them using another parameter A whose 
value is also 1. For the cluster lens the A 1 terms are of numerical order 0.1, A 2 terms are 
of numerical order 0.01 and so on. For the galaxy lens all terms are ~ 1/10 those of the 
cluster. The principal eigenvalues and a r of Eqs. f)A.2j) and (IA.3j) are approximated by 



strong 



A' 



e E D d 



a 



strong 



2 



esc 2 0! <^ cos 3 0i - 2 — esc 2 (f)i cos 0i + 2 Af3 d cos 2 0i + C(A 2 ) 



(20) 



) I esc 2 0i cos 0i - 2 A 2 — esc 4 0i cos 0i + 2 A 3 (cos 0i + (3 d esc 2 0i) + O (A 4 ) \ . 
r d J r d 



These approximate expressions are accurate to at least 0.2% for the strong domain 0.4 0^ < 
0i < 5 0s for the cluster lens and accurate to 0.01% for the galaxy. Strong lensing image 
properties given in equations (ndl -( jl~4"l) differ from conventional Schwarzschild values by 
only a fraction of a percent and are not separately approximated. The effective surface mass 
density k of Eq. f|T2|) . which no longer vanishes as it does for the non-embedded Schwarzschild 
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lens, can be approximated to an accuracy of more than 0.01% for the strong domain as 
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(21) 



An additional strong lensing property of importance is the time delay. It contains a 

-|23|. The arrival time 



AT| P , see 
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geometric part and a potential part, i.e., AT = AT\ 
differences for the two images caused by the difference in geometrical path lengths for our 
embedded Swiss cheese (SC) lens is ATsc| 9 and when computed to maximum accuracy as 



described in 
value AT Sch | 



proves to be almost indistinguishable from the conventional Schwarzschild 



cAT 



SC\g 



cAT 



Sch\g 



[l + z d 




(22) 



The potential part of the embedded lens delay, AT sc | p , as defined in j^, is given by taking 
the difference in the following for the primary and secondary images 



cAT n 



|log 






cot y 



cos 4>\ [ 1 + — cos^ 



d(5 d cos^ 1 + 0{8 2 )\. (23) 



In Figure 6 we have compared AT SC | P = AT p sccondary - AT p primary (i.e., the potential part 
of the time delay) of the cluster lens with the corresponding conventional (non-embedded) 
Schwarzschild value. The reader can see that there is a 2-4% difference in arrival times 
between the theories. 



V. CONCLUSIONS 

This paper is one of a series of investigations of the differences in image properties caused 
by including the gravitational lens's mass in the cosmic mean density. We call such a 
lens an embedded lens. In this paper we have eliminated one of the two impact parameters 
previously required to give embedded point mass lensing quantities such as the bending angle 
a and the lens equation itself. The theory remains more complicated than the conventional 
lensing theory, but is now much easier to use. The new analytical expressions for image 
properties agree with the lowest order results given in jsj]. They can also be compared with 
the higher order results in ij, |2[ that were given as functions of the two impact parameters 
r and 0i. To eliminate r in our prior results for quantities such as a(0i,r o ) in Eq. (32) 
of [l| and obtain results such as Eq. ([7]) given in this paper we had to analytically iterate 
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FIG. 6. A comparison of the potential parts of the time delay (ATsclp) of the embedded cluster 
lens with the conventional theory. The comparison is for sources at same positions even though the 
image positions for the two theories are different. The solid blue line is the conventional potential 
part of the time delay and the dashed red line is the fractional difference of the embedded and the 
conventional theories. 

Eq. (17) of js} to determine ri(^i) and then use the orbit equation (11) of jlj to determine 
r o (0i). The result is given in Eq. flA.lj) of the appendix for completeness and to allow the 
reader to eliminate r in other quantities of interest. 

We have found that with the exception of the potential part of the time delay and the 
effective surface mass density k, strong lensing quantities are only minimally altered by 
making the lens mass a contributor to the mean mass density of the universe. Even there 
the effect is less than 5% on the time-delay for a huge cluster lens, see Fig. 6. For weak lensing 
most effects are also small; however, shear and image ellipticity begin to differ significantly 
(> 10%, see Fig. 5) for large impact angles 0i > 45°. The one quantity that doesn't vanish 
in embedded point mass lensing is k. It turns out to be negative, presumably accounting 
for the missing FLRW mass density in the Kottler void. 

All results given here depend on having a flat (fl = 1) background. Extending them to 
Q 7^ 1 is clearly possible. We expect that many results will differ trivially from what we 
have given here. The applicability of all results given here also depends on the lens being 
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sufficiently condensed so as to be approximated by a point mass. The effects of embedding 



24j. 



on extended lenses remains to be investigated 

To correct for embedding we have used the Swiss cheese cosmologies which are commonly 
criticized for their unrealistic mass distributions, i.e., holes with masses at their centers that 
abruptly appear in otherwise uniform backgrounds. The abrupt discontinuity that appears in 
the cheese is certainly an unrealistic representation of the true matter distribution; however, 
this is primarily an aesthetic complaint. Fortunately for Swiss cheese, its purpose is not to 
represent the mass distribution but instead to account for the effects of mass inhomogeneities 
on the local/global dynamics of the geometry and on the optics of transiting light rays. In 
those two aspects Swiss cheese does quite well. The real shortcoming of a simple Swiss cheese 
type embedded lens (a single condensation moving with the Hubble flow) is the absence of 
any shear at the site of the embedded lens. For such a simple embedded lens, neighboring 
inhomogeneities can only be distributed so as to produce a homogenized gravity field at the 
lens site. Consequently the accuracy of our predictions can be questioned. Stated simply, 
the shortcoming of our lens model, and with standard Swiss cheese itself, is that neighboring 
and distant inhomogeneities produce an homogenized background at the point where the 
lens inhomogeneity is inserted. We suspect this "average" lens is not representative because 
it does not account for effects of local shear. We currently do not have a good estimate 
of how much de-homogenization alters the shielding radius (which is the major source of 
embedding effects) because there are no simple Einstein solutions which accurately model 
local distortions. Such distortions can easily be accommodated in conventional lensing 
theory, but how they would alter the embedding radius is completely unknown. Exact 
Einstein solutions containing a local shear can be constructed by using hierarchical models 
built from Swiss cheese itself. Such a construction will probably be necessary to dependably 
estimate how the spherical shielding radius is distorted and possibly extended by a local 
shear and hence how it modifies predictions made here. 
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Appendix 



The minimum Kottler radial coordinate tq as a function of impact angle 0i (see Fig. 1) is 
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All quantities such as £i, A0, Ar, and p (see Fig. 1 and Fig. 2), previously given as functions 
of 0i and r [l 3 1 can be expressed as functions of the single impact parameter 0i using 

Eq.@3). 

Fig. 2) eigenvalues (a<f,, a r ) 



The azimuthal and radial (with respect to the optical axis, see Fig. 2) eigenval 
of the lensing matrix dO s /d6j to order S 3 as functions of the impact angle 0i and an 
additional lens-geometry parameter (0ED d /r d ) 2 (a term which is of order 6) are 
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